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A general method for the study of the entanglement evolution of graph states under the action of 
Pauli maps was recently proposed in [Cavalcanti et al, Phys. Rev. Lett. 103, 030502 (2009)]. This 
method is based on lower and upper bounds to the entanglement of the entire state as a function only 
of the state of a (typically) considerably-smaller subsystem undergoing an effective noise process 
related to the original map. This provides a huge decrease in the size of the matrices involved in 
the calculation of entanglement in these systems. In the present paper we elaborate on this method 
in detail and generalize it to other natural situations not described by Pauli maps. Specifically, for 
Pauli maps we introduce an explicit formula for the characterization of the resulting effective noise. 
Beyond Pauli maps, we show that the same ideas can be applied to the case of thermal reservoirs 
at arbitrary temperature. In the latter case, we discuss how to optimize the bounds as a function 
of the noise strength. We illustrate our ideas with explicit exemplary results for several different 
graphs and particular decoherence processes. The limitations of the method are also discussed. 

PACS numbers: 03.67.-a, 03.67.Mn, 03.65.Yz 



I. INTRODUCTION 

Graph states [T] constitute an important family of 
genuine niultiparticle-entangled states with several ap- 
phcations in quantum information. The most popu- 
lar example of these are arguably the cluster states, 
which have been identified as a crucial resource for uni- 
versal one-way measurement-based quantum computa- 
tion [21 [3]. Other members of this family were also 
proven to be potential resources very interesting tasks, 
as codewords for quantum error correction [1], to im- 
plement secure quantum communication [S], or to sim- 
ulate some aspects of the entanglement distribution of 
random states [B]. Moreover, graph states encompass 
the celebrated Greenberger-Horne-Zeilinger (GHZ) states 
[7], whose importance ranges from fundamental to ap- 
plied issues. GHZ states can - for large-dimensional sys- 
tems - be considered as simple models of the gedanken 
Schrodinger-cat states, are crucial for quantum commu- 
nication protocols [S], and find applications in quantum 
metrology [9 and high-precision spectroscopy |10j . All 
these reasons explain the great deal of effort made both 
to theoretically understand [T] the properties of, and to 
generate and coherently manipulate, graph states in the 
laboratory [TT] . 

For the same reasons, it is crucial to unravel the dy- 
namics of graph states in realistic scenarios, where the 
system is unavoidably exposed to interactions with its en- 
vironment and/or experimental imperfections. Previous 
studies on the robustness of graph-state entanglement in 
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the presence of decoherence showed that the disentan- 
glement times (i.e. the time for which the state becomes 
separable) increases with the system size [HI [13] . How- 
ever the disentanglement time on its own is known not to 
provide in general a faithful figure of merit of the entan- 
glement robustness: although the disentanglement time 
can grow with the number N of particles, the amount 
of entanglement in a given time can decay exponentially 
with N pL4 . The full dynamical evolution must then be 
monitored to draw any conclusions on the entanglement 
robustness. 

A big obstacle must be overcome in the study of the en- 
tanglement robustness in general mixed states: the direct 
quantification of entanglement involves optimizations re- 
quiring computational resources that increase exponen- 
tially with N. The problem thus becomes in practice 
intractable even for relatively small system sizes, not to 
mention the direct assessment of entanglement during the 
entire noisy dynamics. All in all, some progress has been 
achieved in the latter direction for some very particu- 
lar cases: For arbitrarily-large linear-cluster states under 
collective dephasing, it is possible to calculate the exact 
value of the geometric measure of entanglement through- 
out the evolution [TS]. Besides, bounds to the relative 
entropy and the global robustness of entanglement for 
two-colorable graph states [1] of any size under local de- 
phasing were obtained in Ref. [16) . 

In a conceptually different approach, a framework to 
obtain families of lower and upper bounds to the entan- 
glement evolution of graph - and graph-diagonal - states 
under decoherence was introduced in Ref. [17] . The 
bounds are obtained via a calculation that involves only 
the boundary subsystem, composed of the qubits lying at 
the boundary of the multipartition under scrutiny. This, 
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very often, reduces considerably the size of the matrices 
involved in the calculation of entanglement. No opti- 
mization on the full system's parameter space is required 
throughout. Another remarkable feature of the method is 
that it is not limited to a particular entanglement quanti- 
fier but applies to all convex (bi-or multi-partite) entan- 
glement measures that do not increase under local oper- 
ations and classical communication (LOCC). The latter 
are indeed two rather natural and general requirements 

[Mils]- 

In the case of open-dynamic processes described by 
Pauli maps the lower and upper bounds coincide and the 
method thus allows one to calculate the exact entangle- 
ment of the noisy evolving state. Pauli maps encompass 
popular models of (independent or collective) noise, as 
depolarization, phase flip, bit flip and bit-phase flip er- 
rors, and are defined below. Moreover, one of the vari- 
eties of lower bounds is of extremely simple calculation 
and - despite less tight - depends only on the connec- 
tivity of the graph and not on its total size. The lat- 
ter is a very advantageous property in situations where 
one wishes to assess the resistance of entanglement with 
growing system size. For example, the versatility of the 
formalism has very recently been demonstrated in Ref. 
[20j . where it was applied to demonstrate the robustness 
of thermal bound entanglement in macroscopic many- 
body systems of spin-1/2 particles. 

In the present paper, we elaborate on the details of the 
formalism introduced in |17) . For Pauli maps we give an 
explicit formula for the characterization of the effective 
noise involved in the calculation of the bounds. Further- 
more, we extend the method to the case where each qubit 
is subject to the action of independent thermal baths at 
arbitrary temperature. This is a crucial, realistic type of 
dynamic process that is not described by Pauli maps. In 
all cases, we exhaustively compare the different bounds 
with several concrete examples. Finally, we discuss the 
main advantages and limitations of our method in com- 
parison with other approaches. 

The article is organized as follows: 

• Sec. [nj here we define the notation, introduce def- 
initions and review basic concepts required in the 
following sections. In particular, graph and graph- 
diagonal states are defined in subsection |II A| and 
the noise models considered are presented in sub- 
section III Bl 

• Sec. |III[ a detailed description of the proposed 
framework is given in the context of fully general 
noises. Families of lower and upper bounds for 
the entanglement evolution in the particular multi- 
partition of interest in terms of the entanglement 
of the boundary subsystem alone under an effective 
noise are provided. 

• Sec. |IV| the developed machinery is applied to the 
case of noises described by arbitrary Pauli maps 
and by diffusion and dissipation with independent 



thermal reservoirs at any temperature. Exact re- 
sults for the entanglement decay are obtained for 
Pauli maps, whereas optimized bounds are pro- 
vided in the other cases. 

• Sec. |Vl we first discuss how the method can be 
extended to other initial states or decoherence pro- 
cesses. In particular, how non-tight lower bounds 
for the entanglement evolution of any initial state 
subjected to any decoherence process can be ob- 
tained. Then we comment on the limitations of the 
method. 



• Sec. VI we conclude the paper with a summary of 
the results and their physical implications. 



II. BASIC CONCEPTS 

In this section, we define graph and graph-diagonal 
states, introduce the basics of open-system dynamics and 
the particular noise models used later. 



A. Graph and graph-diagonal states 

Qubit graph states are multiqubit quantum states de- 
fined from mathematical graphs through the rule de- 
scribed below. First, a mathematical graph G(yfi) = 
{V,C} is defined by a set V of iV vertices, or nodes, and 
a set C, of connections, or edges, connecting each node 
i to some other j. An example of such graph is illus- 
trated in Fig. [l] Each vertex i S V represents a qubit in 
the associated physical system, and each edge G C 

represents a unitary maximally-entangling controUed- 
Z {CZ) gate, CZ,, = |0,0j)(0,0,| + |Oaj)(0»lj| + 
|liOj)(liOj| — |lilj)(lilj|, between the qubits i and j 
connected through the corresponding edge. The iV-qubit 
graph state |C?(v,C)o) corresponding to graph G(v,c) is 
then operationally defined as follows: 

(i) Initialize every qubit i in the superposition = 
^(|0i) + |li)), so that the joint state is in the product 

state |g(v)(,) = (g)igv 

(ii) Then, for every connection € C apply the 
gate CZij to |(7(v)g)- That is. 



\G 



(vr)o/ 



C^ijl5(V)o)- 



(1) 



{'J}6C 



Graph state ([T]) can also be defined in an alternative, 
non-operational fashion. Associated to each node i G V 
of a given graph G(v,c) define the operator 



X, 



^0 ' 



(2) 



with Xi and Zj the usual Pauli operators acting respec- 
tively on qubits i and j, and where A/i denotes the set of 
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FIG. 1: (Color online) Mathematical graph associated to a 
given physical graph state. An exemplary bipartition divides 
the system into two subparts: the yellow and white regions. 
The edges in black are the boundary-crossing edges X and 
the nodes (also in black) connected by them are the bound- 
ary nodes y. Together they compose the boundary sub- graph 
G{y,X). The remaining vertices, painted in orange, consti- 
tute the non-boundary subsystem. 



neighbours of i, directly connected to it by an edge {i, j}. 
Operator ([2| possess eigenvalues 1 and —1. It is the i-th 
generator of the stabilizer group and is often called for 
short stabilizer operator. All N stabilizer operators com- 
mute and share therefore a common basis of eigenstates. 
Graph state \G(y c)^) in turn has the peculiarity of being 
the unique common cigenstate of eigenvalue +1 P]. In 
other words, 

^.|G(v,c)o>-|G(v,c)o> V^eV. 

The other 2^ — 1 common eigenstates |G(v,c) ) are in 
turn related to ([T]) by a local unitary operation: 

|G'(v,c)^> = ^^^''\^{Vfi)Q) = Z''\G(yfi)^), (3) 

such that Si \G(y c)ij)—{^'^T^ \G{vfi)v)^ where is a 
multi- index representing the binary string = i^i ... 
with Vi = Q or 1 \/ i € V, and where the short-hand 
notation Z'' = (^jgy Zi"^ has been introduced. There- 
fore, states ([3]) possess all exactly the same entanglement 
properties and, together with |G(v,c)o)) define a complete 
orthonormal basis of H, called the graph-state basis of "H 
(corresponding to the graph G(yfi-)). Any state p diago- 
nal in such basis is called a graph-diagonal state: 

PGD = E^-|G(v,c),)(G(v,c)J, (4) 

where is any probability distribution. Interestingly, 
for any graph, any arbitrary iV-qubit state can always be 



depolarized by some separable map (defined below) into 
the form (|4| without changing its diagonal elements in 
the considered graph basis HJ. 

Two simple identities following from definition ([2| 
will be crucial for our purposes. For every eigenstate 
\G{yfi)^ of S'i, with eigenvalues Si^ = 1 or —1 

UGiy.c),) = (g) ^,|G(v,c)J 

= (g) ^,|G(V.C),>, (5) 

where definition ([2| was used, and 

i"^|G(v,c)J = (-z)Z,.5,® (g) Z,|G(v,c),) 

= s.,(-z)Z, ® (g) Z,|G(v,c),). (6) 

So, when applied to any pure graph ~ or mixed graph- 
diagonal - state, the following operator equivalences hold 
up to a global phase: 

X. O (g Zj, (7a) 
Y,^Z,®^Zj. (7b) 

B. Open-system dynamics 

As we mentioned in Sec. [l] our ultimate goal is to 
study the behavior of graph-state entanglement in real- 
istic dynamic scenarios where the system evolves during 
a time interval t according to a generic physical process, 
which can include decoherence. This process can always 
be represented by a completely-positive trace-preserving 
map A, that maps any initial state p to the evolved one 
after a time t, pt = A(p). In turn, for every such A, there 
always exists a maximum of (B = dim(-H)) operators 
K^j^ such that the map is expressed in the Kraus form 

Operators Ky^ are called the Kraus operators, and de- 
compose the identity operator 1 of in the following 
manner: ^fj_ Kj^K^^ = 1. Conversely, the Kraus repre- 
sentation encapsulates all possible physical dynamics of 
the system. That is, any map expressible as in ([S]) is 
automatically completely-positive and trace-preserving. 
For our case of interest (iV-qubit systems), index p runs 
from to (2^) — 1 = 4^ — 1. For later convenience, we 
will represent it in base 4, decomposing it as the follow- 
ing multi-index: p = pi ... p^, with pi — 0, 1, 2, or 3 

y ieV. 



4 



We call A a separable map with respect to some mul- 
tipartition of the system if each and all of its Kraus 
operators factorize as tensor products of local operators 
each one with support on only one of the subparts. For 
example, if we split the qubits associated to the graph 
shown in Fig. [I] into a set y of boundary qubits (in black) 
and its complement 3^ = of non-boundary qubits (in 
green), A is separable with respect to this partitioning if 
Kfj, = Ky^ (g) Ky^, with Ky ^ and Ky^ operators acting 
non-trivially only on the Hilbert spaces of the boundary 
and non-boundary qubits, respectively. 

In turn, we call A an independent map with respect to 
some multipartition of the system if it can be factorized 
as the composition (tensor product) of individual maps 
acting independently on each subpart. Otherwise, we say 
that A is a collective map. Examples of Jully independent 
maps are those in which each qubit i is independently 
subject to its own local noise channel £i. By the term 
independent map without explicit mention to any respec- 
tive multipartition we will refer throughout to fully inde- 
pendent maps. In this case, the global map A factorizes 
completely: 

K{p)^£i®£2® ...^SNip). (9) 

It is important to notice that all independent maps are 
necessarily separable but a general separable map does 
not need to be factorable as in ^ and can therefore be 
both, cither individual or collective. 

1, Pauli maps 

A crucial family of fully separable maps is that of 
the Pauli maps, for which every Kraus operator is pro- 
portional to a product of individual Pauli and iden- 
tity operators acting on each qubit. That is, = 

... Miv) t^lMi'X'- ■ •«'CriV^„ = ^/P^J.(^|M, with (Tio = 

(the identity operator on qubit i), an = Xi, — Yi, 
and crj3 = Zi, and P(^i, ... = any probability 
distribution. Popular instances are the (collective or 
independent) depolarization and dephasing (also called 
phase damping, or phase-flip) maps, and the (individ- 
ual) bit-flip and bit-phase-flip channels (22j. For exam- 
ple, the independent depolarizing (D) channel describes 
the situation in which the qubit remains untouched with 
probability 1 — p, or is depolarized - meaning that it is 
taken to the maximally mixed state (white noise) - with 
probability p. It is characterized by the fully-factorable 
probability = pi^^ x ... pn^,^, with p^^ = I - p 
and pii — Pi2 = Pi:} — p/3, V i G V. The indepen- 
dent phase damping (PD) channel in turn induces the 
complete loss of quantum coherence with probability p, 
but without any energy (population) exchange. It is also 
given by a fully-factorable probability with piQ = l—p/2, 
Pii = = pi2, and = p/2, V i G V. 

For later convenience, we finally recall that each Pauli 
operator Ui^. can be written in the following way: 



Ti = Z^\X'^' , with Ui and = 0, or 1. Indeed, 

notice that <Ji2vi+\vi+ui\.2 = Ti(«i.t)i) (up to an irrelevant 
phase factor for m = 1 — Vi), where "| I2" stands for 
modulo 2. In this representation, called the chord rep- 
resentation [.23 , the Kraus decomposition of the above- 
considered general Pauli map has the following Kraus op- 
erators: Kc(u,v) = V^ci 

TN{uM,vr^) = ^/Pc(u.y) T(uy), where U = (wi, ... un) 
and V = (ui, ... vn). The probability Pc{u.v) = 
Pc(ui,vi, ... Miv,ujv) in turn is related to the original P^ by 

Pc(ui,Vi, ... Miv.Uiv) = P{2vi + \vi+Ui\2, ... ,2VK-i-\vN-i-UN\2)}- 

2. Independent thermal baths 

An important example of a non-Pauli, indepen- 
dent map is the generalized amplitude-damping channel 
(GAD) 122] ■ It represents energy diffusion and dissipa- 
tion with a thermal bath into which each qubit is indi- 
vidually immersed. Its Kraus representation is 



^.M,=2 ^ ^^_r^(v/r^|0.)(0.| + |1.)(1,|), (10c) 
and 

^.p.=3^y^=^|l,:>(0.|- (lOd) 

Here n is the average number of quanta in the thermal 
bath, p = p{t) = 1 — e~2T(2"+i)* -g ^j^g probabihty of the 
qubit exchanging a quantum with the bath after a time 
t, and 7 is the zero-temperature dissipation rate. Chan- 
nel GAD is actually the extension to finite temperature 
of the purely dissipative amplitude damping (AD) chan- 
nel, which is obtainen from GAD in the zero-temperature 
limit n = 0. In the opposite extreme, the purely diffu- 
sive case is obtained from GAD in the composite limit 
n 00, 7 — > 0, and 717 = F, where F is the diffusion 
constant. Note that in the purely-diffusive limit, channel 
GAD becomes a Pauli channel, with defining individual 
probabilities p^g = 1(1 - p/2 + y/l -p), = 2 ^ p^^, 

andp,3 = l(l-p/2- VT^), V z G V. 

Finally, the probability p in channels D, PD and GAD 
above can be interpreted as a convenient parametrization 
of time, where p — Q refers to the initial time and p = \ 
refers to the asymptotic limit < — > 00. 

III. EVOLUTION OF GRAPH-STATE 
ENTANGLEMENT UNDER GENERIC NOISE 

As mentioned before, the direct calculation of the en- 
tanglement in arbitrary mixed states is a task exponen- 
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tially hard in the system's size In this section, we 
elaborate in detail a formalism that dramatically sim- 
plifies this task for graph - or graph-diagonal - states 
undergoing a noisy evolution in a fully general context. 
Along the way, we also describe carefully which require- 
ments an arbitrary noisy map has to satisfy so that the 
formalism can be applied. 



A. The general idea 

Consider a system initially in graph state ([T]) that 
evolves during a time t according to the general map 
(Isl) towards the evolved state 



Pt 



A(|G(v,c)o» = E^A< 



|G(v,c)o)(G'(v,c)ol^i- 



(11) 



where \g(y)g) = ®i^y l+i)- other words, we explicitly 
factor all the CZ gates corresponding to non-boundary 
qubits out. 

Consider now the application of some Kraus op- 



erator of a general map on graph state (12 1 

\g(^\^)- The latter can al 



ways be written as {^^^ j}e7? 
with 



K„ 



cz,,J,,y^i, (13) 



Now, consider every map A such that transformation rule 
( 13 1 yields, for each fi, modified Kraus operators of the 



form 



We would like to follow the entanglement E{pt) of pt 
during its entire evolution. Here, E is any convex en- 
tanglement monotone [TSl [TU] that quantifies the entan- 
glement content in some given multi-partition of the sys- 
tem. An example of such multi-partition is displayed in 
Fig. [T] where the associated graph is split into two sub- 
sets, painted respectively in yellow and white in the fig- 
ure. The edges that connect vertices at different subsets 
are called the boundary-crossing edges and are painted in 
black in the figure. We call the set of all the boundary- 
crossing edges X CC, and its complement X = C/X the 
set of all non-boundary-crossing edges. All the qubits 
associated to vertices connected by any edge in X con- 
stitute the set 3^ C V of boundary qubits (or boundary 
subsystem), and its complement y = V /y is the non- 
boundary qubit set. We refer to G(^y^x) a-s the boundary 
sub-graph. 

We can use this classification and the operational def- 
inition ([I]) to write the initial graph state as 



\G 



(V,C) 



CZ,j\G(y^X)Q/ 



K. 



(14) 



where Ky^ and Ky^ are normalized modified Kraus op- 
erators acting non-trivially only on the boundary and 
non-boundary qubits, respectively. In the last, 7 = 
{pi, i & y} and uo — {pi, i G y} are multi-indices la- 
beling respectively the alternatives for the boundary and 
non-boundary subsystems, being pi in turn the individ- 
ual base-4 indices introduced after Eq. Q . The modified 
map A, composed of Kraus operators is then clearly 
bi-separable with respect to the bi-partition "boundary 
/ non-boundary". For all such maps the calculation of 
E{pt) can be drastically simplified, as we see in what 
follows. 



In these cases, the evolved state (11) can be writen as 



Pt = A(|G(v.c)o)) 



GZi. 



Pt Q$ CZki. 

{k,i}ex 



(15) 



)®\9iy)^)^ (12) with 



{i,3}ex 



Pt = l{\G(y^x)^) ® l5(3J)o)) ^Y.^y-^M\^(y-^)^)^'^iy-^)^\^i(i.)®^y^{p)\9m 

w 7 

I 



(16) 



where Ky(~^\^-^ is the 7-th modified Kraus operator on the 
boundary subsystem given that Ky^ has been applied to 
the non-boundary one. Recall that both 7 = 7(/i) and 
Lo = a;(/x) come from the same single multi- index p and 
are therefore in general not independent on one another. 
In the second equality of ( 16 ) we have chosen to treat w 



as an independent variable for the summation and make 
7 explicitly depend on u. This can always be done and 
will be convenient for our purposes. 

The crucial observation now is that the GZ operators 
explicitly factored out in the evolved state (15 1 corre- 
spond to non-boundary-crossing edges. Thus, they act 
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as local unitary operations with respect to the multi- 
partition of interest. For this reason, and since local uni- 
tary operations do not change the entanglement content 
of any state, the equivalence 



E{pt) = E{pt) 



(17) 



holds. 

In the forthcoming subsections we see how, by exploit- 
ing this equivalence in different noise scenarios, the com- 
putational effort required for a reliable estimation (and 
in some cases, an exact calculation) of E[pt) can be con- 
siderably reduced. The main idea behind this reduction 
lies on the fact that, whereas in general state (111 the 



entanglement can be distributed among all particles in 
the graph, in state ( 16 1 the boundary and non-boundary 
subsystems are explicitly in a separable state. All the 
entanglement in the multi-partition of interest is there- 
fore localized exclusively in the boundary subgraph. The 
situation is graphically represented in Fig. [2j where the 
same graph as in Fig. [l] is plotted but with all its non- 
boundary-crossing edges erased. 





More precisely, the general approach consists of ob- 
taining lower and upper bounds on E{pt) by bounding 
the entanglement of state (16) from above and below as 
explained in what follows. 



1. Lower hounds to the entanglement evolution 

The property of LOCC monotonicity of E, which 
means that the average entanglement cannot grow during 

an LOCC process [TH], allows us to derive lower-bounds 
on E{pt). The ones we consider can be obtained by the 
following generic procedure: 



after bringing the studied state into the form 
( |16[ ), apply some local general measurement M = 
with measurement elements M^^, on the 
non-boundary subsystem 3^; 



(ii) for each measurement outcome oj trace out the mea- 
sured non-boundary subsystem; 



(Hi) and finally, calculate the mean entanglement in the 
resulting state of the boundary subsystem 3^, aver- 
aged over all outcomes uj. 



FIG. 2: (Color online) Same graph as in Fig. [Tjbut where all 
non-boundary-crossing edges have been erased, representing 
the fact that the boundary and non-boundary subsystems are 
fully unentangled. The entanglement in the whole system is 
obtained via a calculation involving only the smaller boundary 
subsystem. 



Since this procedure constitutes an LOCC with respect 
to the multipartition under scrutiny, the latter average 
entanglement can only be smaller than, or equal to, that 
of the initial state, i.e. : 



E{pt) = E{p,) >Y.P^E[Y.-{g^y^^\ky^M^^^^^^^ (18) 

UJ uj' 7 

I 

with = X;^,(5(5;)^|7^1^,Mt.M<^i^5;^,|.g(5j)^) being the Notice that if the states {Ky^,\g(y)^)} of 
probability of outcome u. the non-boundary subsystem happen to be or- 
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thogonal, then there exists 
surement M = {M^ = 



an 



optimal mea- 



that 
that 



can distinguish them 



, l3(y)o)<9(y)ol^y^ 

unambiguously, so 

= 5, 

t 



and P,., = {g^y^jK'-^KyJg^y^^ 



} 



X 



In these cases an optimal lower bound is achieved as 



E{pt) > 

! Ky{y\uj)\G^y^X)^) (G(3;,A')o ■ (19) 



Full distinguishability of the states in the non-boundary 
subsystem allows to reduce the mixing in the remaining 
boundary subsystem. In other words, the measurement 
outcome w works as a perfect flag that marks which sub- 
ensemble of states of the boundary subsystem, from all 
those present in mixture (16 1, corresponds indeed to the 
obtained outcome. 

In the opposite extreme, when states {-l^yu'ldQj)^)} 
are all equal, no flagging information can be obtained 
via any measurement. In this case, the resulting bound 
is always equal to that obtained had we not made any 
measurement at all, but just directly taken the partial 
trace over y from (16): 



E{pt) > 



\G 



'{y,x}o){G{y,x)„ 



k 



(20) 



where |3^| stands for the number of non-boundary qubits 
and full mixing over variable lo takes place now. 



Henceforth we refer to lower bound ( 20 ) as the lowest 



lower hound (LLB). As its name suggests, its tightness 
is far from the optimal one given by (19). However, as 



we will see in the forthcoming subsections, due to the 
partial tracing, it typically does not depend on the total 
system's size but just on the boundary subsystem's. 

This constitutes an appealing, useful property, for it 
allows one to draw generic conclusions about the robust- 
ness of entanglement in certain partitions of graph states, 
irrespective of their number of constituent particles (see 
examples in Sec. IV B). 



2. Upper bounds to the entanglement evolution 



On the other hand, we consider upper-bounds on E{pt) 
based on the property of convexity of E, which essentially 
means that the entanglement of the convex sum is lower 
than, or equal to, the convex sum of the entanglements 
[mils]. From (fiel), the latter implies that 



E{pt) = E{pt)<Y,P^E 



1 

P 



(21) 



where, once again, P^ = '^9(y)^}K^y^kyJg(y)^). In 
each term of the last summation the boundary and 
non-boundary subsystems inside the brackets are in 
a product state. Therefore, as for what the multi- 
partition of interest concerns, the non-boudary sub- 
system works as a locally-added ancila (in a state 
'h^y^\9(y)^){9(y)^\K^yJ and consequently does not 
have any influence on the amount of entanglement. This 
leads to the generic upper bound 



3. Exact entanglement 



Notice that upper bound ( 22 ) and optimal lower bound 



(19) coincide. This means that, in the above-mentioned 
case when states {Ky^\g(y-^ )} are orthogonal, these co- 
incident bounds yield actually the exact value of E{pt): 



K 



E{pt) 

t 



y{i\'^)\G{y,x)o){G(y,x)„\Kyf^^^^^). (23) 



P.,E{ 



K 



E{pt) < 

y(i\^)\G{y,x)o){G{y.x)„\Kl^„^^^.^). (22) 



Expression (23) is still not an analytic closed formula 



for the exact entanglement of pt , but reduces its calcula- 
tion to that of the average entanglement over an ensem- 
ble of states of the boundary subsystem alone. More in 
detail, a brute-force calculation of E(pt) would require 
in general a convex optimization over the entire (2^) - 
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complex-parameter space. Through Eq. ( 23 ) in turn 



that condition ( 14 1 is trivially satisfied. Since for every 



such calculation is reduced to that of the average en- 
tanglement over a sample of 2l-^l states (one for each u) 
of |3^| qubits, being |3^| the number of boundary qubits. 
The latter involves at most 2'-^' independent optimiza- 
tions over a (2l-^l) -complex-parameter space. This, from 
the point of view of computational memory required, ac- 

2 

counts for a reduction of resources by a factor of (2l^l) . 
Alternatively, when computational memory is not a ma- 
jor restriction - for example if large classical-computer 
clusters are at hand -, one can take advantage of the 
fact that the |3^| required optimizations in (23) are in- 



dependent and therefore the calculation comes readily 
perfectly-suited for parallel computing. In this case, 
it is in the required computational time where an |3^|- 
exponentially large speed-up is gained. 

In the cases where states {KyJ\g(y-^^} are not orthog- 
onal and the upper and lower bounds do not coincide. 



expressions ( 22 1 and ( 18 ) still yield highly non-trivial up- 



per and lower bounds, respectively, as we discuss in Sec. 

ttvH 

Finally, it is important to stress that all the bounds de- 
rived here are general in the sense that they hold for any 
function fulfilling the fundamental properties of convex- 
ity and monotonicity under LOCC processes. This class 
includes genuine multipartite entanglement measures, as 
well as several quantities designed to quantify the useful- 
ness of quantum states in the fulfillment of some given 
task for quantum-information processing or communica- 
tion ^9j. 



IV. GRAPH STATES UNDER PAULI MAPS OR 
THERMAL RESERVOIRS 

In the present section we apply the ideas of the pre- 
vious section to some important concrete examples of 
noise processes. This shows how the method is helpful 
in the entanglement calculation for systems in natural, 
dynamic physical scenarios. We first discuss the case of 
Pauli maps and then the generalized amplitude damp- 
ing channel (thermal reservoir). Explicit calculations for 
noisy graph states composed of up to fourteen qubits are 
presented as examples. 



A. Pauli maps on graph states 



Pauli maps defined in Sec. II B provide the most im- 
portant and general subfamily of noise types for which 
expression ( 23 1 for the exact entanglement of the evolved 



state applies. In this case, every Xi or Yi Pauli matrix 
in the map's Kraus operators is systematically substi- 
tuted by products of Zi and 1^ according to rules Q. 
The resulting map A defined in this way automatically 
commutes with any CZ gate and is fully separable, so 



qubit in the system four orthogonal single qubit opera- 
tors are mapped into products of just two, several differ- 
ent Kraus operators of the original map contribute to the 
same Kraus operator of the modified one. This allows us 
to simplify the notation going from indices /ij , which run 
over 4 possible values each, to modified indices fli hav- 
ing only two different alternatives. In fact, the original 
operators give rise to only 2^ modified ones of the 
form 



(24) 



where multi-index /x stands for the binary string /2 = 
jji ... p,N, with /Ltj = or 1, Vi G V. Probability is 
given simply by the summation of all in the original 
Pauli map over all the different events fi for which cr^ 
yields - via rules ([t]) - the same modified operator Z'^ in 



(24) 



To compute the latter modified probability we move 
to the chord notation [23 , mentioned at the end of Sec. 
II Bl Indeed, under transformation ([t]), we have that 



5<8>jeM^J%SO thatr(i;,y) 



l(ui,t)i)'< 



Pa 



The latter coincides with Z^ every time [li = \vi 



V i G V. Thus, in this representation, the 



modified probability Pcp. is obtained from the defining 
probability Pccuy) in ttis original map by the explicit 
formula 



P 



U 



,"jv,IAi-IZjs 



(25) 

The modified Kraus operators (p4l) in turn are fully 



separable; thus, as said, they trivially satisfy factoriza- 
tion condition (14). We can express them as = 
with 



= \l P(^\Cj)Z~' and 



'PiZ". 
P6) 

The new multi- indices are 7 = {/Xi, i S y} and uj — 
{/ij, * S y}, and the corresponding probabilities satisfy 



-^(7|lj)-R'J 




states {Ky 




^63^75(|0») + (-l)''1l■')) = ^/^-'l9(3^)-,)} are 
trivially checked to be all orthogonal. Thus, they provide 
perfect flags that mark each sub-ensemble in the bound- 
ary subsystem's ensemble. The perfect flags are revealed 
by local measurements on the non-boundary qubits in 
the product basis Therefore, for Pauli maps 

the exact entanglement E{pt) can be calculated by ex- 
pression (23), which, in terms of binary indeces 7 and 
UJ, and using graph-state relationship ([3|, can be finally 
expressed as 

Eipt)^J2P^^(Y.Pi7\^)\^(y.xu){G(^y,x)J), (27) 
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14 qubits 



turn: 




0.1 0.2 0.3 0.4 0.5 0.6 



FIG. 3: (Color online) Negativity vs p for 14-qubit (green tri- 
angles) and 12-qubit (pink solid curve) cluster states under- 
going independent depolarizing noise, and for the bipartition 
shown in inset. Parameter p can be thought as a parametriza- 
tion of time (see text). 



B. Independent thermal reservoirs on graph states 

In the case of Pauli maps the entanglement lower and 
upper bounds coincide, and deliver the exact entangle- 
ment. However, this is not the case for general, non- 
Pauli, noise channels. The upper bound is given, as 
usual, by convexity. The lower bounds must be optimized 
by appropriately choosing the LOCC operations. Here, 
we investigate and optimize measurement strategies for 
channel GAD, defined in Sec. [ITB21 



Observe that the Kraus operators defined in Eqs. ( 10 ) 
satisfy the following: and commute with any CZ 
operator, while for every j S V different from i it holds 

'ij.{Kii (g) Zj) and {Ki^ (g) 
Based on this, one can 



that {Kii (g) lj).CZij 
'^j)-CZij = CZij.(K, 



i3 



^ CZ, 
g) Z,) 



perform the factorization in equation il3h and apply this 



way the formalism described in Sec. HI 



mX 



In what follows we focus on two main limits of channel 
GAD discussed in Sec. |IIB 2| the purely-dissipative limit 
n = (amplitude damping) , and the purely-difusive limit 
n — ^ oo, 7 — >■ 0, and n7 = T. 



1. Graph states under zero-temperature dissipation 



In Fig. [3] we have plotted the bipartite entanglement of 
the exemplary bipartition of one qubit versus the rest 
shown in its inset for fourteen and twelve qubit graph 
states evolving under individual depolarization. This 
map, as said before, is characterized by the one-qubit 
Kraus operators ^/l — pi, y^p/3X, ^Jj^Y , and yJp/ZZ. 
The parameter p (0 < p < 1) refers to the probabil- 
ity that the map has acted: for p = the state is left 
untouched and for p = 1 it is completely depolarized. 
Once more, p can be also set as a parametrization of 
time: p = Q referring to the initial time (when nothing 
has occurred) and p = 1 referring to the asymptotic time 
t — > 00 (when the system reaches its final steady state). 

As the quantifier of entanglement, we choose the nega- 
tivity [25] , defined as the absolute value of the sum of the 
negative eigenvalues of the density matrix partially trans- 
posed with respect to the considered bipartition. It is a 
convex entanglement monotone that in general fails to 
quantify entanglement of some entangled states - those 
ones with positive partial transposition (PPT) - in di- 
mensions higher than six |19j . However, since its cal- 
culation does not require optimizations but just matrix 
diagonalizations, it is very well-suited for a simple illus- 
tration of our ideas. 

We emphasize that, for the graph used in Fig. [3j 
a brute-force calculation would involve diagonalizing a 
X 2^'' = 16384 X 16384 density matrix for each value 



We consider a four qubit linear (ID) cluster state sub- 
jected to the AD map and study the decay of entangle- 
ment in the partition consisting of the first qubit versus 
the rest shown in the inset of Fig. |4j Along with the 
exact calculation of entanglement via brute-force diago- 




of p, whereas with the assistance of expression (27 1 E(pp) 
is calculated via diagonalization of many 2^ x 2r — 8x8 
dimensional matrices only. 



0.2 0.4 



FIG. 4: (Color online) Negativity vs p for a 4 qubit lin- 
ear cluster subjected to the amplitude damping channel and 
for the partition displayed in the inset. Solid (upper) black 
curve: exact entanglement; Solid (lower) blue curve: lowest 
lower bound LLB (obtained by tracing out the flags) ; Dashed- 
dotted blue curve: upper bound (obtained by convexity); Dot- 
ted red curve: LB(0) (obtained by measuring the flags in the 
Z basis); Dashed brown curve: LB(7r/4) (obtained by mea- 
suring the flags in the X basis). 
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FIG. 5: (Color online) Lower bound LB{9) to the negativity 
as a function of the angle 9 in the measurement basis, for the 
same situation as in Fig. |4] and for fixed values of p. Each 
value p = 0.01,0.1,0.3,0.5 has two curves associated to it. 
The horizontal (gray) straight line represents the exact en- 
tanglement at each p, while the blue (black) curve represents 
the bound LB(6') at this p. The red line (vertical) corresponds 



to 61 = f , i.e. 



measurements in the basis {]+), ]— )}■ 



nalization of the partially-transposed matrices, the lowest 
lower bound LLB (20), obtained by tracing out the flags, 



and the upper bound ( 22 1 , obtained from convexity, are 
plotted. In addition, the tightness of the lower bounds 



(181 obtained by the flag measurements can be scanned 



as a function of the measurement bases. 

Based on observations about the behavior of the sys- 
tem under the AD map we can guess good measurement 
strategies. For example, examination of the initial state 
reveals that at p = each of the non-boundary qubits 
is in one of the states of the basis {|+), |— )}; whereas at 
p = 1, in one of the states of {|0), |1)}. We call the lower 
bound obtained from ( |18[ ) through measurements in the 
basis {|+), |— )} LB(7r/4), and the one obtained from ([l8| 
through measurements in {|0), |1)} LB(0). The latter 
bounds are the two additional curves plotted in Fig. |4] 
We observe that LB(0) provides only a slight improve- 
ment over the LLB, whereas LB(7r/4) appears to give a 
significant one. This raises the obvious question of how to 
optimize the choice of measurement basis at each instant 
p in the evolution. 

As an illustration we consider lower bounds LB(0) 
obtained from ( |18[ ) through orthogonal measurements 
composed of projectors j^-l-) — cos6'|0) + sin^jl) and 
16*—) = — sin^jO) -|-cos0|l), and look for the angle 6 that 
gives us approximately the largest value of LB{9). This is 
certainly not the most general measurement scenario one 
may consider, but it gives one a hint on how to increase 
the tightness of the bounds. 



FIG. 6: (Color online) Same as Fig. |5]for p = 0.9. 



Figs. [5] and [6] illustrate this idea. At fixed values of p, 
we have varied angle 9 in the range [0, 7r/2]. The entan- 
glement given by LB{9) for each value of 9 is compared 
with the exact entanglement at the given p. In phys- 
ical terms, we are taking snapshots of the evolution of 
the system's entanglement at discrete time instants. The 
value of 9 at each instant p that maximizes LB(^) repre- 
sents the optimal measurement basis at that particular 
instant. As is clearly seen in Fig. [5] for small values of p 




0.2 0.4 



FIG. 7: (Color online) Negativity vs p for a 4-qubit linear 
cluster for the bipartition shown in the inset, subjected to 
the generalized amplitude damping channel in the diffusive 
limit n — >■ oo. The central curve corresponds to both the ex- 
act entanglement (solid black) and LB(7r/4) (dashed brown), 
which coincide exactly. The lower curve corresponds to both 
LBB (solid blue) and LB(0) (dotted red), which also coincide 
exactly. The upper curve is the upper bound ( 22 1 (dot-dashed 
blue). 
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angles around 9 — j give the closest approximations to 
the exact entanglement, in consistence with the signifi- 
cant improvement of LB(7r/4) over the LLB observed in 
Fig. |4j For large values of p though, the best approxi- 
mations tend to be given by the angles away from 6 — 
as can be observed in Fig. [6j It must still be kept in 
mind that none of these closest approximations equals 
the exact entanglement of the state. 



2. Graph states under infinite-temperature difusion 

We now consider the purely-diffusive case of the GAD 
channel, where each qubit is in contact with an indepen- 
dent bath of infinite temperature. In Fig. [7] we display 
the entanglement evolution in a similar way as in Fig. [4] 
Since in the purely-diffusive limit channel GAD becomes 
a Pauli map, as was mentioned in the end of Sec. II B 2[ 
bound LB(7r/4) yields the exact entanglement. LB(0) 
on the other hand coincides with the lowest lower bound 
LLB. The fact that in this case LB (6') reaches the exact 
entanglement at — j can also be seen in a clearer way 
in Fig. m 

In Fig. [t] upper bound ( 22 1 is plotted as well. Since in 
this case the channel is a Pauli channel, one would expect 
the upper bound to coincide with the exact entanglement 
as well. The fact that this does not occur is because, 
even though the noise itself is describable as a Pauli map. 



> 
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FIG. 8: (Color online) Lower bound LB{9) to the negativity 
as a function of the angle 9 in the measurement basis, for 
the same situation as in Fig. [Tj and for fixed values of p. 
Each value p — 0.01,0.1,0.3,0.5 has two curves associated to 
it. Again, the horizontal (gray) straight hne represents the 
exact entanglement at each p, while the blue (black) curve 
represents the bounds LB(0) at the this p. The red line (ver- 



tical) corresponds to 
{|+),|-)}- 



i.e. measurement in the basis 



the plotted upper bound has been calculated using the 
original Kraus decomposition of Eqs. (10 1, which is not 
in a Pauli-map form. For every given particular Kraus 
decomposition of a superoperator, the naive application 
of convexity always yields UB through Eq. (22 1, but this 



needs not the tightest, for the Kraus decomposition of a 
superoperator is in general not unique. This observation 
leads to a whole family of upper bounds for a given map. 
In the same spirit as with the lower bounds, one could 
in principle optimize the obtained UBs over all possible 
Kraus representations of the map. 



V. EXTENTIONS AND LIMITATIONS 

The framework developed here is not restricted to 
graph states. The crucial ingredient in our formalism 
are to factor out all the entangling operations that act 
as local unitary transformations with respect to the con- 
sidered partition, and to redefine of the Kraus operators 
acting on the state, reducing the entanglement evaluation 
problem to the boundary sub-system alone. Given an 
entangled state and a prescription for its construction in 
terms of entangling operations, useful bounds and exact 
expressions for the entanglement can be readily obtained. 
As an example, a GHZ-like state = a |0)^^-|-/3 |1)®^ 
can be operationally constructed by the sequential ap- 
plication of maximally-entangling operation CNOTij = 

(l,:Oj| to 



|0,0j) (O.Ojl + |0,1,) (0,1,1 |liO,-) (lil.l -f |1,1, 



the product state (a| 



/3|1)) ® |0) (g) ... ® |0) such that 



IV-) = (g)CA^Or,,,+i(a|0) + m) ® |0) ® ... ® |0). Us- 

ing our techniques and the permutation symmetry of the 
state it can be seen that, for GHZ-like states as above un- 
dergoing the previously discussed noise processes, the en- 
tanglement evaluation in any bipartition can be reduced 
to that of a two qubits system. It is also important to 
mention that the techniques presented here can also be 
straightforwardly extended to higher-dimensional graph 
states |27. 

In addition, it is important to mention that all bounds 
developed so far can in fact also be exploited to follow the 
entanglement evolution when the system's initial state is 
a mixed graph-diagonal state. This is simply due to the 
fact that any graph-diagonal state as ^ can be thought 
of as a Pauli map Aqjj acting on a pure graph state: 

PGD = ^P,y\G(Vfi)J{G{V,C)J = 

V 

^P.znG(v,c)o)(G(v,c)jZ''-AGi5(|G(v,c)o))- (28) 

Thus, the entanglement at any time i in a system initially 
in a mixed graph-diagonal state pgd; and evolving under 
some map A, is equivalent to that of an initial pure graph 



state \G 



whose evolution is ruled by the composed 



map A o Agd, where h.Qu is defined in (28 1. When A 
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is itself a Pauli map, then A o Kgd is also a Pauli map 
and the expression (27) for the exact entanglement can 



be applied. For the cases where A is not a Pauli map but 



the relations ( 13 ) are satisfied by its Kraus operators, the 



relations ( 13 ) will also be satisfied by the composed map 



A o Kgd, so that all other bounds derived here also hold. 

Furthermore, as briefly mentioned before, any arbi- 
trary state can be depolarized by some separable map 
towards a graph-diagonal state without changing the di- 
agonal elements in the considered graph basis [21]. The 
latter, since the entanglement of almost all states cannot 
increase under separable maps |28j . implies that all the 
lower bounds presented here also provide lower bounds 
to the decay of the entanglement that, though in general 
far from tight, apply to almost any arbitrary initial state 
subject to any decoherence process. 

The gain in computational effort provided by the ma- 
chinery presented here diminishes with the ratio between 
the number of particles in the boundary subsystem and 
the total number of particles. For example, for multipar- 
titions such that the boundary system is the total system 
itself, or for entanglement quantifiers that do not refer to 
any multi-partition at all, our method yields no gain. An 
example of the latter are the entanglement measures that 
treat all parties in a system indistinguishably, some of 
which, as was mentioned in the introduction, have been 
studied in Refs. [TSIIIS]. These methods naturally com- 
plement with ours to offer a rather general and versatile 
toolbox for the study of the open-system dynamics of 
graph-state entanglement. 



VI. CONCLUSIONS 

We have studied in detail a general method for comput- 
ing the entanglement of graph and graph-diagonal states 



undergoing decoherence, introduced in Ref. [17] . This 
method allows to drastically reduce the effort to com- 
pute the entanglement evolution of graph states in sev- 
eral physical scenarios. We have given an explicit formula 
for the construction of the effective noise involved in the 
calculation for Pauli maps and extended the formalism 
to the case of independent baths at arbitrary tempera- 
ture. Also, we have elaborated the formalism to construct 
non-trivial lower and upper bounds to the entanglement 
decay where exact results cannot be obtained from the 
formalism itself. 



Finally we would like to add that the necessary require- 
ments on the noise channels for the method to apply do 
not prevent us from obtaining general results for a wide 
variety of realistic decoherence processes. Furthermore, 
the conditions required on the entanglement measures 
are satisfied by most quantifiers. 
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